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We report on the scaling of transport properties around the vortex melting in YBa2Cu307_,s 
oriented-twin single crystals in applied magnetic fields between IT and 18T. We find that for all the 
measured field range the linear resistivity scales as p(t,0) ~ t sy J r ±(sin(9)t~ sx ), with t = \T — Tbc\ 
and the angle between de planar defects and the magnetic field. The scaling is valid only for 
angles where the transition temperature Tbg(@) shows a cusp. The critical exponents sx and sy are 
in agreement with the values predicted by Lidmar and Wallin only at magnetic fields below 4T. A 
change in the value of sx from sx = 1 ± 0.2 to sx — 3 ± 0.2 at around H cr ~ 4T when the magnetic 
field is increased, is responsible for changes in the shape of the Tbg(0) curve and in the dependence 
of the linear dissipation on temperature and angle. The results strongly suggest the existence of a 
different vortex glassy phase in twinned crystals compared to the Bose-glass state found in samples 
with linear defects. 



One of the most important aspects in high tempera- 
ture superconductors is the existence of a melting line in 
the magnetic field-temperature (H — T) vortex phase di- 
agram. The quest to understand the thermodynamic na- 
ture of the transition has been reflected in an important 
number of publications!!! Nowadays it is accepted that 
in clean crystals the solid-liquid transition is first order, 
while the introduction of random point-like or correlated 
disorder transforms it into second order. Due to the high 
value of the Ginsburg parameter in high T c materials, 
the critical region associated with this kind of transition 
becomes experimentally accessible. This important fact 
has led to the development of scaling theoriesBB which de- 
scribe the vortex dynamics near the transition, enclosing 
its properties into two critical exponents. In both the- 
ories the critical exponents settle the way in which, the 
correlation length, £ ~ \T — Tbg\~ v — t~ v , and the time 
scale in which a fluctuation of this size survives, r ~ £ 2 , 
diverge as the transition is approached. In the Bose-glass 
theory linear correlated defects induce an anisotropic cor- 
relation volumefl with £|| ~ £^ which are, respectively, 
the correlation lengths parallel and perpendicular to the 
correlated defects. In this theoryB the linear resistivity is 
predicted to scale as 

p(t,Q) ~t ay F ± (sm(0)t- ax ) (1) 

with a phase boundary separating the Bose-glass and the 
vortex liquid of the form 

T BG (0) = T BG (0)[1 - (5E^)V~] ( 2 ) 

Xc 

where sy = v[z — 2), sx = 3w, and x c = sin(0 c )t sx . 
Lidmar and Wallin pointed outa that the correct ther- 
modynamic variable to be included in the scaling anzats 
for the dynamical response of the vortex system ua tilted 
magnetic fields is B, instead of H as used in refu. This 



important remark leads to a modification of eqs. |lj and g 
through a change in the dependence of sx on the critical 
exponent vnp Accordingly, the form of the new phase 
boundary is given by eq. @ with sx — v, while the linear 
resistivity still scales as eq. |l| but again with sx — v. 

Equation || describes one of the key characteristics of 
the Bose-glass to liquid transition, i.e. the existence 
of a cusp in the transition curve at = Q. This cusp 
was first observed by Worthington kL alu in twinned 
YBa 2 Cu307_(5 crystals. Fleshier et alH studied the vor- 
tex dynamics for different orientations of the twin bound- 
aries with respect to the applied field and current, and 
explained the cusp in terms of pinning effects by the 
twin potentials. The relation between the cusb and the 
Bose-glass transition was shown by Jiang et ala through 
transport measurements in YBa2Cu307_5 crystals with 
columnar defects. Later on, Grigera et ala identified 
a Bose-glass phase in twinned YBa2Cu30y_5 crystals 
through a scaling analysis of the vortex dynamics near 
the cusp line. An interesting point related to the shape 
of the cusp, however, has gone unnoticed in the litera- 
ture. While in samples with columnar defects the tran- 
sition curve js— spJl described by Tbg(0) ~ 1 — \0\ at 
small anglesil I'i'J 11 ! in twinned crystals a non-linear tran- 
sition curveBQ&HLij is observed. Monte Carlo simula- 
tions of the Bose-glass meltingu in a vortex system in- 
teracting with linear defects gave v ra 1 which, through 
eq. ||, seems to account only for the results in samples 
with columnar defects. It is therefore important to clar- 
ify whether the cusp in samples with planar correlated 
defects arises as a consequence of a liquid to Bose-glass 
transition as suggested in refE, and to which extent the 
scaling theory describing the vortex dynamics is applica- 
ble to this more anisotropic case. 

In this paper we show that within the linear regime 
the scaling hypothesis p(t,6) ~ t sv J r ±(sin(8)t~ sx ) is ful- 
filled for all magnetic fields between IT and 18T in 
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YBa2Cu307-«5 heavily oriented-twin single crystals. For 
all applied fields the transition temperature Tbg shows 
a cusp around 6 — 0. However the shape of the Tbg(&) 
curve is only consistent with the predictions of Lidmar 
and Wallin for H < H cr ss 4T. Above this field a change 
of the cusp's shape is accompanied by a sudden increase 
in sx from sx = 1 ± O.f to sx = 3 ± 0.2 while sy re- 
mains constant in the whole field range. This jump of sx 
is clearly reflected in the scaling of the linear resistivity. 
We discuss these findings in terms of a change in the uni- 
versality class of the vortex solid-liquid transition which 
we adscribe to a corresponding change in the nature of 
the vortex solid phase as the field is swept through H cr . 

The YBa2Cu307_,5 crystal with oriented twins was 
mounted onto a rotatable sample holder with an angu- 
lar resolution of A9 « 0.05° inside a cryostat provided 
with an 18T magnet. The crystal has a critical temper- 
ature of T c ~ 93K and AT C ~ 0.5K width dimension of 
1 x 0.5 x 0.03mm 3 . A uniform dc-current was injected 
through the crystal parallel to the twin planes and the 
crystal was rotated at angles 9 in a plane perpendicu- 
lar to the twins in a constant Lorenzt force configuration 
(see inset Fig. |l|). In this setup the vortices are moved 
and tilted in a direction perpendicular to the defects. All 
transport measurements were performed within the lin- 
ear regime at J=3Acm~ 2 . 

In Fig. [I] we show the transition temperature Tbg($) 
for different magnetic fields. This temperature has been 
defined from p(9, T) measurements as the temperature 
at which the data falls below our experimental resolu- 
tion (V/I < 20/ifi). The value at 9 = coincide within 
experimental error with that obtained from the scaling of 
the linear resistivity. It is seen that the shape of the cusp 
changes as the magnetic field is increased. The change is 
rather abrupt, and occurs at a crossover field H cr w AT. 
Below H cr the Tbg{8) is linear in \9\, while above a non- 
linear dependence develops. The solid lines, which are fits 
to eq. H with sx = 1 for H= IT and 2T, and sx = 3 for 
6T and 8T illustrate this point. The arrows correspond 
to a critical angle, 9 Cl above which the scaling properties 
of the linear resistivity are lost. Not only the shape of the 
cusp changes, but its amplitude increases with field. The 
angular width of the cusp shrinks when the field increases 
but remains almost constant for H > 4T. Interestingly 
enough, the cusp exists up to the highest magnetic fields 
we measured (18T). This may be compared to the behav- 
ior repotted for Tbg{0) m (Ba,K)BiC>3 irradiated cubic 
crystals,E9 where the cusp gradually disappears upon sur- 
passing the matching field, B^, in which the number of 
vortices equals the number of the columnar defects. The 
absence of saturation in our experiments indicates that 
the twin boundary potentials are still efficient to induce 
vortex confinement and localization in a field range in 
which induced columnar tracks with « 18T would 
destroy the superconducting properties of the crystal. 

The changes in the shape of the phase boundary are 
also reflected in the dissipation in the liquid state as 
shown in Fig. g, in which p(T) is plotted for different 



angles near 9 = 0. Panel (a) shows data characteristic 
of the low field regime, while panel (b) displays the data 
corresponding to the high field region. It is worth notic- 
ing that the linear dissipation has a different behavior at 
low and high fields as the angle between the magnetic 
field and the defects is increased. At low H the p(T) 
curves shift almost parallel to each other while at high 
fields p(T) stretches to lower temperatures as 9 is in- 
creased. These differences are indicative of a change in 
the vortex dynamics below and above a crossover field 
H cr which separates the low and high field regions. In 
order to quantify the observed changes we refer to the 
scaling theoryEm according to which the linear resistiv- 
ity, when the applied field is parallel and vortices move 
perpendicular to the correlated disorder, should scale as 
p{T,0) ~ t sy with sy = v(z - 2). On the other hand, 
the resistivity as a function of the tilting angle at the 
transition temperature should behave as p(0, 9) ~ 9 sy / sx 
for small angles. These scalings are shown in the insets 
of panels (a) and (b) for the corresponding fields. In 
general we find that for H < H cr , sy = 2.4 ± 0.3 and 
sy/sx = 2.3 ± 0.2, while for H > H cr sy = 2.7 ± 0.3 and 
sy/sx = 0.9 ± 0.2. Note that within experimental error 
sy seems to be field independent while the ratio sy/sx 
undergoes an abrupt change. From these we conclude 
that the changes in the behavior of the linear resistivity 
are solely due to a change in the exponent sx. In order 
to see this, we show in Fig. || the scaled linear resistivity 
as predicted by eq. [j] Panel (a) shows the same data dis- 
played in Fig. ||(a), while in panel (b) the scaling of the 
data in Fig. ||(b) is shown. Both scalings are very good 
for p(T, 9) curves below the corresponding critical angle 
9 C . Curves for 9 > 9 C do not scale. An important point 
to note here is that the critical exponent sx is different at 
low and high fields. In panel (a), sx = 1 ± 0.2 was used 
to obtain a good scaling, while in panel (b) the value of 
sx was increased to 3 ± 0.2 in order to make the curves to 
scale. In the insets the same scaling is shown but forcing 
the value of sx to be 3 at low fields and 1 at high fields, 
and maintaining the rest of the parameters constant, as 
required by consistency with eq. [j] at 9 = 0. Clearly, this 
assignment of critical exponent values renders the scaling 
impracticable in both cases, indicating that the change in 
the value sx is outside the errors involved in such scaling. 
It is worth noticing that an exponent sx — 1 corresponds, 
through eq. 0, to a linear Tbg{®) dependence, which in- 
deed is seen in Fig. [l] for low fields, while a value of 
sx = 3 gives a non-linear shape as observed in Fig. ^ at 
high fields. 

In Fig. ^ we show a summary of our findings by plot- 
ting the values of the two critical exponents and their 
ratio as a function of the magnetic field. The open cir- 
cles in panel (a) correspond to the values of sy obtained 
by a fit of the resistivity data to p(t, 0) ~ t sy , from which 
also Tbg(0) was obtained. With these values of Tbg(0) 
the scaling of p(t, 9), as required by eq. |], was performed 
to obtain the solid squares in panel (a) and (b). The 
error bars were estimated by the effect the variation of 
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the exponents sx and sy has on the scaling of the linear 
resistivity. As is seen sy is, within error bars, indepen- 
dent of the magnetic field and of the way in which it is 
obtained. The value of sx, however, jumps from 1 to 3 
at 4T, remaining constant for higher or lower fields. Fi- 
nally in panel (c) we show the ratio of sy/sx. The solid 
squares were calculated from the values in panel (a) and 
(b) while the open circles were obtained by adjusting the 
resistivity data at T BG (0) to p(Q,9) ~ 9 sy l sx . Note that 
the change at around 4T is independent of the way in 
which sy/sx was obtained. 

An important conclusion that can be drawn from our 
results is that the way in which the linear dissipation 
in the vortex liquid varies as a function of angle and 
temperature can be accounted for by a scaling law of the 
form of eq. [l] independently of the applied magnetic field. 
This is a fingerprint of a second order vortex solid-liquid 
phase transition occurring at Tbg{9) for angles below 
C . Within this physical picture, our results indicate the 
existence of two different field regimes characterized each 
one by different values of the critical exponents which 
govern the shape of the transition curve and the dynamic 
response of the liquid phase in the linear regime. On 
general grounds the two regimes can be interpreted in 
terms of a change in the universality class to which the 
solid-liquid transition belongs. Let us first discuss the 
low field regime. 

Our results for H < H cr closely resemble those ob- 
served in crystals with columnar defects. We may there- 
fore assume that the Tbq(9) is linear in \0\ as a con- 
sequence of the interaction of the vortex system with 
some kind of columnar-type potentials. The identifica- 
tion of the columnar-type defects which might be present 
in our twinned crystal is however difficult. A possi- 
ble candidate we found by optical inspection of our 
crystal were end points of twin planes located between 
the voltage contacts. An end point breaks the transla- 
tional symmetry of the twin plane and may emulate a 
columnar- type defect for vortices located near it. This 
end point corresponds microscopically to zones where a 
twin turns into a transverse one, or more likely in a mi- 
crotwinned region which can not be seen under an op- 
tical microscope. Other source of columnar-type defects 
which can also contribute to localise the vortices in this 
regime, might be screw dislocations which also exist in 
twinned YBa2Cu307_a crystals.E.3 Within this assump- 
tion a matching field should exist, above which all lin- 
ear defects become occupied by a vortex line. In view 
of the work of Klein et al. showing that the cusp in 
Tbg(9) exists up to magnetic fields more than three times 
larger than the corresponding B^, it is not straightfor- 
ward in our case to identify the crossover field H cr with 
the matching field of the linear defects. It could well be 
that H cr represents the field above which the effects of 
the columnar defects are completely washed out. In this 
case H cr should be sample dependent and could even be 
absent in crystals of very high quality. This physical pic- 
ture could be tested in a more direct way by measuring 



the vortex transport properties in an oriented-twin crys- 
tal with a low concentration of columnar tracks induced 
by ion irradiation in the c-axis direction. 

Once the effect of the columnar potentials has been sat- 
urated, the dynamics of the vortex system is determined 
by its interaction with the twin planes. This corresponds 
to the field regime above H cr . In this field range vor- 
tices localise inside the twin potentials, at Tbg(0)- How- 
ever, the vortex solid that is formed above H cr ought to 
be different from that present in samples with columnar 
tracks alone. A way to see this is to consider the form 
in which vortices can accommodate inside the correlated 
potentials. While columnar tracks allow single vortex 
occupancy only, in a finite sample with twin planes vor- 
tices will continue to accommodate inside the correlated 
defects until the increase in energy due to their mutual 
interaction overcomes the energy gain due to the vortex- 
defect interaction. This difference in the way in which 
vortices interact with and become localised at the twin 
potentials gives rise to a different dynamical response of 
the vortex liquid compared to that found in a system 
interacting with columnar tracks. It is therefore reason- 
able to conclude that the vortex liquid will freeze in a 
vortex solid which does not have the same structure as 
a Bose-glass. The non-linear phase boundary in ob- 
served in our experiments may be taken as a fingerprint 
of this novel vortex solid phase. In general we may say 
that the form of the phase boundary separating the solid 
and liquid phases is determined by the type of correlated 
disorder which is present in the sample. 

In summary, we have shown that the vortex solid to 
liquid transition in YBa2Cu30y_5 oriented-twin crystals 
is second order within IT and 18T, and occurs at a phase 
boundary which is non-linear in angle when the vortex 
dynamics is dominated by the interaction with planar de- 
fects. The linear dissipation in the liquid state in tilted 
magnetic fields can be very well described by a scaling law 
theoretically proposed to explain the dynamics of a vor- 
tex system interacting with columnar-type pinning po- 
tentials, but with different critical exponents. These ob- 
servations strongly indicate the existence of a distinctly 
different vortex solid phase in twinned crystals compared 
to the Bose-glass phase found in samples with columnar 
defects. 
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FIG. 1. Angular dependence of Tbg(9) - T BG (0). Solid 
lines are fits to eq. ^| with sx = 1 for data at IT and 2T, and 
sx = 3 for data at 6 and 8T. The arrows indicate the critical 
angle 9 C above which the linear resistivity does not scale. Data 
point are mirrored, gray points correspond to mirrored data. 
Inset: experimental setup for transport measurements. 
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Fig. 2: B, Moiorovand E. Osquiguil 

FIG. 2. (a) Normalized resistivity p — p/p(95K) for dif- 
ferent tilting angles up to 9 C at H=1T, 9 =0.2, 0.45, 0.65, 
1.0, 1.3, -1.3, 1.6, 1.9, 2.2, 2.7, 3.2, 3.7, -4.5, 6.0,-7.0, 8.0, 
-8.0 deg. Upper inset: p(t, 0) . The solid line is a fit using 
p ~ t sy with T SG (0) = (89.82 ± 0.05)K and sy = 2.5 ± 0.3. 
Lower inset: ,5(0,0). The solid line is a fit to p ~ 9 ay/3X 
with sy/sx = 2.2 ± 0.3. (b) idem for H=9T, 9 =-0.15, 0.3, 
0.75, 0.85, 1.1, 1.25, 1.47, -1.8, 2.25 degrees. Upper inset: 
Tbg(0) = (79.85 ± 0.1)K and sy = 2.8 ± 0.3. Lower inset: 
sy/sx = 1.1 ±0.3. 
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FIG. 3. (a) Scaling of the resistivity curves shown in 
figure ^(a) according to eq. |l| for sy = 2.3, sx — 1.1 and 
Tbg = 89.82K. The inset shows the same scaling using sx = 3. 
(b) Scaling of the resistivity curves shown in figure ^](b) ac- 
cording to eq. [j] for sy = 2.8, sx = 3.0 and Tbg = 79.85K. 
The inset shows the same scaling using sx — 1. 
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FIG. 4. Summary of the results of the scalings for different 
fields, (a) Values of sy according to the resistivity scaling 
(solid squares) and obtained from p(t,0) ~ t sy (open circles), 
(b) Values of sx according to the resistivity scaling, (c) val- 
ues of sy/sx from the fitting of /5(0, 8) (open circles) and the 
ratio of sy and sx obtained from the solid squares in panels 
(a) and (b). 
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